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Abstract. We present results of a study of the soft X-ray excess 
ijtijthe black hole candidate Cygnus X-1 made with the Rosat 
^SPC, using observations taken during persistent emission at 
®bital phases close to 0.5. The soft excess can be well fitted as a 
^ckbody with temperature kTbb = 0.13 ± 0.02 keV. kTbb did 
©tjvary appreciably with intensity of the source. By assuming 
fljat the distance of the source is its lower limit of 2.5 kpc, a lu- 
Tfflnosity of the soft excess of 4.7 • 10^'^ erg was obtained. 
*^^m this, disc temperatures were calculated as a function of 
radius, assuming the compact object to be a 10 Mq black hole, 
^^articular, the temperature at 7 Schwarzschild radii expected 
wbe highly representative of the total disc emission. This was 
and to be 0. 1 3 keV, in very good agreement with the spectral 
ling result. This good agreement strongly supports the identi- 
l^njation of the soft excess with emission from the disc around a 
®ack hole. 
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1. Introduction 

Cygnus X-1 is a high mass X-ray binary consisting of a 09.7 
supergiant and a compact object which is one of the best can- 
didates for a black hole, based on determinations of the mass 
function and lack of X-ray eclipses (Paczyfiski 1974). Two pe- 
riodicities are seen in the data: the 5.6 d orbital period known 
from optical spectroscopy (Gies & Bolton 1982) and a 294 d 
periodicity in the X-ray light curve thought to be related to pre- 
cession of the accretion disc (Priedhorsky et al. 1983). Dips in 
the X-ray intensity have been observed in the source at orbital 



Send offprint requests to: M. Balucinska-Church 



phases close to zero, thought to be due to absorption in inhomo- 
geneities in the stellar wind from the companion (Kitamoto et 
al. 1989; Balucinska & Hasinger 1991). 

The X-ray spectrum is complex, consisting of a number of 
components. The high energy spectrum shows a break at ~ 200 
keV and has been described by Comptonisation models (eg 
McConnell et al. 1994). A strong spectral steepening below 10 
keV corresponding to transitions to the high X-ray state of the 
source has been found (Tananbaum et al. 1972). 

The presence of a reflection component was identified in 
Exosat GSPC data by Done et al. (1992) who found that an ab- 
sorption edge at ~ 6 keV and a broad bump in the spectrum 
peaking at ~ 20 keV could be well fitted as a reflection term 
similar to those seen in AGN. This is a very broad spectral term 
which contributes to the continuum between 0.1 and ~ 100 
keV. Additionally there is an iron emission line at ~ 6.2 keV 
(Barret al. 1985). 

At lower energies the observed spectrum has a soft X- 
ray excess which became apparent through inconsistencies be- 
tween the column density Nh obtained from the low energy cut- 
off of the X-ray spectra and the galactic Nh ~ 6 ± 2 • 10^^ H 
atom cm~^ derived from the reddening towards the compan- 
ion HDE 226868. These inconsistencies were removed if the 
X-ray spectrum is modelled by a two-component model con- 
sisting of a blackbody to represent the soft excess plus a power 
law to represent the Comptonisation term in the band 1-20 keV. 
The soft excess was first detected by Priedhorsky et al. (1979). 
Based on the Einstein/MPC data, Balucifiska and Chlebowski 
(1988) showed that this soft excess may contribute up to 5% 
of the flux in the 1.2-20 keV band. A steep low energy spec- 
trum and an emission feature at 0.73 keV were found in Exosat 
transmission grating spectrometer (TGS) data in the band 0.4 
- 2.0 keV by Barr and van der Woerd (1990). The presence of 
the soft excess was confirmed by extensive studies of the Ex- 
osat ME observations by Balucinska and Hasinger (1991), who 
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Fig. 1. X-ray light curve for the Rosat observation at phase 0.46 in the energy band 0. 1 - 2.0 keV with 10 s timebins. The dotted horizontal lines 
correspond to the intensity ranges in which spectra were accumulated 



showed that the excess could be modelled as a blackbody with 
temperature kTbb ~ 0.25 keV. However, the ME data were not 
able to constrain the temperature very well, as the peak of the 
blackbody spectrum was outside the range of the instrument. 

The Rosat PSPC allows a proper study of the soft excess 
in the appropriate energy band containing the peak of the emis- 
sion. We present detailed spectral fitting analysis of the soft ex- 
cess in the band 0.1- 2.0 keV, from observations made when 
Cygnus X-1 was in its low state. Results are given here for the 
part of the observation at orbital phase close to 0.5 where there 
are no dips and at which the effect of extra absorption due to the 
stellar wind of the companion is expected to be minimal. Re- 
sults from the rest of the observation in which intensity dips are 
seen will be given in a second paper. 

2. Data analysis 

The observations were made on 1991 April 18-20 spanning a 
total of 66 h and consisting of three sections of continuous ob- 
servation lasting 2030 s, 1890 s and 1100 s. The first two sec- 
tions were close to orbital phase zero and we do not discuss 
this data here. The third section was at orbital phase 0.46 calcu- 
lated using the accurate ephemeris of Gies and Bolton (1982). 
Source data were extracted from a circle of radius 0.2 1 ° centred 
on the source to include photons scattered outwards in the dust 
halo (Predehl & Schmitt 1995). Background was extracted from 
an annulus extending from 0.22° to 0.29° . The background was 
subtracted and the data were corrected for deadtime and vi- 
gnetting. We present results below obtained using the March 
1992 response function appropriate to this observation in the 
early part of the Rosat mission (see Fiore et al. 1994). However 
in order to assess systematic errors, we also analysed the data 
using the January 1993 response, and found that parameters de- 
rived from spectral fitting changed by only 1.5% to 3%, giving 
confidence in the results. 



The light curve of the third section of data is shown in Fig. 
1. It can be seen that, because of the strong variability of the 
source on short timescales, it is not possible to accumulate data 
for a sufficiently long time with the source at a constant inten- 
sity level to produce a spectrum of suffficient quality (requiring 
~ 1 00 s of data) . Consequently the data were divided into seven 
intensity bins, each 10 c/s wide covering the ranges from 60 - 
70 c/s to 120 - 130 c/s in the band 0.1 - 2.0 keV. These intensity 
bins were equivalent typically to 130 s of data, so that Poisson 
errors were sufficiently small for good fitting to be made. For 
each of these intensity bands, a spectrum was produced by re- 
binning into 24 channels. We added conservatively a system- 
atic error of 2% to each channel. If the actual error is somewhat 
smaller, the effect will be that the errors we give for parameters 
derived from spectral fitting should be reduced. The values of 
reduced would have to be increased by 10 - 20%, so that un- 
acceptable models become more unacceptable; however for our 
best fit models x'i would remain less than unity. 

Additionally, a total spectrum was produced including all 
intensities with 24 channels to allow searching for emission and 
absorption features. 

3. Results 

Firstly, spectral fits of the spectra in the 7 intensity bins were 
tried using simple models including an absorbed power law 
model. In this case, the quality of the fits was unacceptable, 
with reduced typically as large as 1.24 with 21 degrees of 
freedom. Moreover, the power law photon index T averaged 
over the 7 intensity bands was 2.3 ± 0.2 (Icr), and the individ- 
ual values of T in all 7 bands were significantly higher than 
the widest range of possible values we obtained in the higher 
Exosat ME energy band, between 1.5 and 1.8 (Baiucifiska & 
Hasinger 1991). This clearly demonstrates the presence of the 
soft excess. The Rosat PSPC is not able to constrain power law 



the band 1-10 keV or higher is needed. It is however very well 
able to determine the parameters of a thermal component peak- 
ing within the PSPC band, and so the most sensible approach is 
to use the PSPC data in conjunction with T values previously 
established for the source. For this we have used results from 
observations of Cygnus X-1 with the TTM and HEXE instru- 
ments onboard Mir. These instruments are normally used simul- 
taneously and so provide high quality spectra extending over 
the very broad 1 - 200 keV band. Consequently continuum pa- 
rameters may be well established and we have previously found 
that from several TTM/HEXE observations T for the intrinsic 
power law has the value 1.73 ± 0.02 (Balucinska-Church et al. 
1996). 

Fitting the simple absorbed power law model with T fixed 
at 1 .73 gave very poor fits in all 7 intensity bands with Xr typ- 
ically 1.7. Adding a second soft power law to the model to pa- 
rameterise the spectra was next tried, and in all cases this re- 
quired r ~ 8, clearly unphysical. The large inconsistency be- 
tween this value and 1 .73 is further evidence for the soft excess. 
Thus we were able to discard a simple one-component model, 
and adopted a physically more realistic two-component model 
consisting of a power law with T fixed at 1 .73, plus a blackbody 
to model the soft excess, with a single absorption term for both 
components. Good fits were obtained in all 7 intensity bands, 
with reduced significantly improved over simple models. 



isation clearly increasing linearly with intensity. 

The mean value of kTbb for all intensity bands was 
. 1 3 ± . 2 ke V, showing that the PSPC was able to locate the 
blackbody peak and so determine the value of kTbb with good 
accuracy. It was not necessary to specify a separate column 
density for each component, and the mean value of Nh was 
7.5 ± 1.0 • 10^1 H atom cm-2. 
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Fig. 2. (a) Spectral fitting results for the 7 intensity levels into which 
the observation was divided showing blackbody temperature kTbb as a 
function of c/s. The dashed hue shows the mean value, (b) Variation of 



the power law normalisation at 1 keV Ipi (photons cm 



keV"^) 



and the blackbody normaUsation Lbb (expressed as luminosity at 10 
kpc in erg s~^) with intensity 



In Fig. 2a, results for the blackbody temperature kTbb as 
a function of count rate are shown. It can be seen that there is 
no variation of kTbb with intensity, showing that the short-term 
variability apparently takes place without change in the shape 
of the spectrum. Only the normalisations of the power law and 
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Fig. 3. Top panel: Spectral fit to the data in the total spectrum using 
the two-component model, together with residuals. Bottom panel: in- 
dividual model component contributions to the unabsorbed total en- 
ergy spectrum. The dotted and dashed lines show the power law and 
blackbody components respectively 



We tested the robustness of our result for kTbb by fitting the 
total spectrum, adding all intensity bands, with the same two- 
component model, and fixing F at values covering the widest 
range of acceptable values: 1.5 to 1.8. For these extreme val- 
ues of F, kTbb changed from 0.15 ± 0.02 keV to 0.13 ± 0.02 
keV; ie not significantly different from the value we previously 
obtained in the individual intensity bins. 

In Fig. 3 we show the best fit to the total spectrum together 
with residuals with photon index F fixed at 1.73; the Individ- 



in the lower panel. We see no evidence for emission or absorp- 
tion features, particularly the feature at ~ 0.7 keV reported by 
Barr and van der Woerd (1990). Best fit values from fitting the 
total spectrum are: Nh = 7.0 ± 1.0 • 10^^ H atom cm~^, black- 
body bolometric luminosity at 10 kpc: L = 7.5lg 3 • 10^^ erg 
s~^, kTbb = 0.14tg keV, power law normalisation at 1 keV 
= 1.37 ± 0.12 photons cm'^ s"! keV^ (90% confidence er- 
rors). It can be seen that the blackbody luminosity is non-zero 
at a high confidence level; the value is greater than 2 • 10^^ erg 
at 99% confidence. The absorbed flux of the blackbody in 
the band 0.1 - 2.0 keV is 2.3 • 10~^° erg s~^, the power law ab- 
sorbed flux in the same band is 9.0 • 10~^° erg s~^, so that the 
blackbody comprises 20% of the energy flux in this band. For 
unabsorbed spectra, this fraction rises to about 50%. 

4. Discussion 

We have investigated further the lack of any spectral feature at 
~ 0.7 keV, reported by Barr and van der Woerd (1990). Using 
their best fit results we have simulated the Exosat TGS data, and 
find that if the feature was present at the same strength in our 
Rosat observation, it would have been easily detectable. How- 
ever the PSPC data had a 25% higher flux in the band below 0.7 
keV, but a 25% smaller flux above 0.7 keV; ie there were signif- 
icant changes in both the softness and brightness of the source. 

The results consistently show a temperature for the soft ex- 
cess of ~ 0.13 keV. We now estimate the temperature of the 
emission expected from an accretion disc about a 10 Mq black 
hole, and show that there is good agreement. 

It is expected that the disc thermal emission will be modified 
by reflection; however the work of Ross and Fabian (1993) has 
shown that for reflection in the discarounda IOMq black hole, 
and in the restricted band 0.3 to 3 keV (Fig. 5, Ross & Fabian) 
the incident blackbody is not strongly modified. Furthermore, 
Ross et al. (1992) showed that the temperature of the accre- 
tion disc at 7 Schwarzschild radii, where the contribution to the 
luminosity r^F peaks (F is the emitted flux), is highly repre- 
sentative of the total disc emission from the multi-temperature 
disc. They also showed that fitting a synthesised disc spectrum 
with a simple blackbody produces results for kTbb in reason- 
able agreement with the temperature at 7rs . 

Firstly we calculate the luminosity of the soft component 
at the time of observation, using our best fit spectral param- 
eters and setting the distance of the source at its lower limit 
of 2.5 kpc (Bregman et al. 1973). This was found to be Lx = 
4.7 • 10^'^ erg s~^. However the actual distance of the source 
is poorly known and the luminosity could be several times 
larger. This luminosity was converted into mass accretion rate 
M ~ 5.2 • 10^'^ g , assuming an efficiency of 0.1, and 
from this we calculated disc temperatures as a function of radius 
in a multi-temperature disc (eg Frank, King & Raine 1992). At 
7 Schwarzschild radii (210 km), we find a temperature kTbb of 
. 1 3 ± . 3 ke V, where the errors reflect only the uncertainties 
in the soft component luminosity from spectral fitting. Thus, 
there is good agreement of the observationally determined soft 



tance of the source is appreciably greater than 2.5 kpc, the lumi- 
nosity will be higher, however temperatures in the disc are pro- 
portional to M^/^ and so a factor of 2 in distance implies only 
a 40% increase in temperature. 

The blackbody luminosity is less than 1% of the Edddington 
luminosity Le = 1.47 • 10^^ erg s"^ for a 10 Mq object with 
70% hydrogen. The luminosity of the power law component is 
estimated as ~ 1 .9 • 10^^ erg s~ ^ Thus the above comparison 
of the temperature of the soft excess with disc temperatures is 
for the low state of the source at a low fraction of the Eddington 
limit. 

Using the total spectrum, the emitting area of the blackbody 
was found to be 1.3 •10^'^ cm^ from the luminosity of this com- 
ponent and the temperature kTbb =0.14 keV. If this is emission 
from a disc, the outer radius would have to be ~ 460 km, so that 
the emission region effectively extends from the inner disc to ~ 
15 Schwarzschild radii. 
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